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We investigate a new type of approximation to quantum 
determinants, the "quantum Fredholm determinant" , and test 
numerically the conjecture that for Axiom A hyperbolic flows 
such determinants have a larger domain of analyticity and 
better convergence than the Gutzwiller-Voros zeta functions 
derived from the Gutzwiller trace formula. The conjecture is 
supported by numerical investigations of the 3-disk repeller, 
a normal-form model of a flow, and a model 2-d map. 
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I. INTRODUCTION 

While various periodic orbit formulas may be for- 
mally equivalent, in practice some are vastly preferable 
to others. Trace formulas, such as the thermodynamic 
averages in classical dynamics, and the semi-classical 
Gutzwiller trace formula in quantum mechanics are dif- 
ficult to use for anything other than the leading eigen- 
value estimates. However, dynamical zeta functions [jlj, 
Fredholm determinants Q and the Gutzwiller-Voros zeta 
function ||E| have recently been established as powerful 
tools for evaluation of classical and quantum averages 
in low dimensional chaotic dynamical system s |5| - |9|. 
Recent advances include new, cycle expansion JlOj] based 
numerical spectra evaluations and Riemann conjecture 
inspired functional equations. The convergence of cy- 
cle expansions of zeta functions and Fredholm determi- 
nants depends on their analytic properties; particularly 
strong results exist for Axiom A hyperbolic systems, for 
which the dynamical zeta functions are meromorphic [ pT| , 
and the classical Fredholm determinants are entire func- 
tions |l2|jr^ ]. No such results exist for the determinants 
used in quantum theory, but formal analogies to the clas- 
sical case have led to introduction of the quantum Fred- 
holm determinant 
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as an alternative to the Gutzwiller-Voros zeta function 
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We present here the numerical evidence in support of the 
conjecture of ref. 14 1: 



*PC: version 16 july 1993, file chaos.nbi.dk: 
/ users /predrag/tex / articles / q_fredholm / qfred6 .tex 



For Axiom A systems the quantum Fredholm determinant 
has a larger domain of analyticity than the Gutzwiller- 
Voros zeta function. 

We shall consider here only purely hyperbolic flows 
with the topology of a Smale horseshoe. The important 
conceptual insight of Smale |l5j is the realization that for 
such flows the associated zeta functions have nice analytic 
structure. In a more formal setting, such flows are called 
"Axiom A" , and Ruelle [0 proves that for expanding an- 
alytic maps the zeta functions are meromorphic, and the 
spectrum is discrete. This differs very much from the 
intuition acquired in studies of quantum chaos; there is 
no "abscissa of absolute convergence" and no "entropy 
wall" , the exponential proliferation of orbits can be con- 
trolled, and the Selberg-type zeta functions are entire and 
converge everywhere. 

In classical mechanics and number theory, the zeta 
functions and the Fredholm determinants are exact. All 
quantum mechanical studies take a saddle-point approx- 
imation - the Gutzwiller trace formula - as the starting 
point, and for quantum mechanics an important con- 
ceptual problem arises already at the level of deriva- 
tion of the Gutzwiller-Voros zeta functions; how ac- 
curate are they, and can the periodic orbit theory of 
such semi-classical approximations be systematically im- 
proved? We shall not address this problem here; in this 
paper we are interested in the convergence of cycle expan- 
sion truncations of the Gutzwiller-Voros zeta function. 



1 



The problem is classical in the sense that all quantities 
used in periodic orbit calculations - actions, stabilities, 
Maslov phases - are classical quantities. 

The main limitation of the study presented here is that 
the Fredholm determinants are proven to be entire only 
for hyperbolic flows with symbolic dynamics with finite 
grammar (Axiom A flows). Hence our numerics is re- 
stricted to three systems with complete binary grammar 
and bounded nonlinearity: 

(a) the 3-disk repeller 

(b) a Hamiltonian Henon map as a normal form ap- 
proximation to a flow 

(c) a model 2-d hyperbolic mapping 

The paper is organized as follows: in sect. |n| we review 
the e volu tion operator formalism for smooth flows. In 
sect. Ill we explain the theorems that guarantee that 
Fredholm determinants for Axiom A systems are entire. 
In sect. IV we define and motivate the determinants used 
in this paper. In sect. ^ we review the cycle expansions, 
and give convergence estimates for various expansions. 
In sect. VI we discuss the numerical evidence in support 



A. Trace formula for maps 

If the evolution is given by a discrete time mapping, 
and all periodic points are known to have stability eigen- 
values Afc ^ 1 strictly bounded away from unity, the 
trace C n is given by the sum over all periodic points x of 
period n: 
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of the quantum Fredholm determinant conjecture. 



is the [dxd] Jacobian matrix evaluated at the periodic 
point x, and the product goes over all periodic points Xi 
belonging to a given prime cycle p. The trace formula is 
the Laplace transform of tr£* which, for discrete flows, 
is simply the generating function 



II. FLOWS, EVOLUTION OPERATORS AND 
THEIR SPECTRA 

Functional determinants and zeta functions arise in 
classical and quantum mechanics because in both the 
dynamical evolution can be described by action of linear 
evolution operators on infinite-dimensional vector spaces. 
The classical evolution operator for a d-dimensional map 
or a (d + l)-dimcnsional flow is given by: 



C t (y,x) = S(y-f(x))g t (x) 



(3) 



For discrete time, f n (x) is n-th iterate of the map /; 
for continuous flows, f t (x) is the trajectory of the initial 
point x. g t {x) is a weight multiplicative along the tra- 
jectory; its precise functional form depends on the dy- 
namical average under study. For purposes of this sec- 
tion it suffices to take g t (x) — 1, essentially the Perron- 
Frobenius operator case. 

The global averages (escape rates, energy eigenvalues, 
resonances, fractal dimensions, etc.) can be extracted 
from the eigenvalues of the evolution operators. The 
eigenvalues are given by the zeros of appropriate deter- 
minants. One way to evaluate determinants is to expand 
them in terms of traces, log det = tr log, and in this way 
the spectrum of an evolution operator becomes related to 
its traces, it. periodic orbits. Formally, the traces tr£* 
are easily evaluated as integrals of Dirac delta functions 
as follows: 
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where e "°, e Vl , e " 2 , . . . are the eigenvalues of C. For 
large times det (l — jW(ij)) — > Aj, where Aj is the 
product of the expanding eigenvalues of so the 

trace is dominated by 
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and diverges at the leading eigenvalue 1/z — e~ v ° . This 
approximation, which in current physics literature is 
called the "thermodynamic" or the "/ of a" formal- 
ism |l6| , is adequate (but far from optimal) for extraction 
of the leading eigenvalue of C, and difficult to apply to 
extraction of the non-leading eigenvalues. 



B. Trace formula for flows 

For flows the eigenvalue corresponding to the eigenvec- 
tor along the flow (the velocity vector) necessarily equals 
unity for all periodic orbits, and therefore the integral 
([|) requires a more careful treatment . 

To evaluate the contribution of a prime periodic orbit 
p of period T p , one choses a local coordinate system with 
a longitudinal coordinate dxu along the direction of the 
flow, and d transverse coordinates x± 



dx ± dx {l S(x ± - fi(x))6( X]] - /f(a:)) ■ (7) 
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Integration is restricted to an infinitesimally thin tube V p 
enveloping the cycle p. 

Let v = \F(x)\ be the velocity along the orbit, and 
change the longitudinal variable to dx || —vdr. Whenever 
the time t is a multiple of the cycle period T p , the integral 
along the trajectory yields 

~ oo „ 

/ dx\\5(x\\ - /|f(aO) = ~ rT v) / dT 

Jv„ r=1 Jp 



6(t - rT p ) 
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Linearization of the periodic flow in a plane perpendic- 
ular to the orbit yields the same weight as for the maps 

1 



/ dxj_5(xj_ 
Jv„ 



fI rTp (x)) 



dct 1 



(9) 



where J p is the p-cycle [dxd] transverse Jacobian, and 
we have assumed hyperbolicity, ie. that all transverse 
eigenvalues are bounded away from unity. A geometri- 
cal interpretation of weights such as @ is that after the 
r-th return to a surface of section, the initial tube V p 
has been stretched out along the expanding eigendirec- 
tions, with the overlap with the initial volume given by 
l/|det(l- J;) L 

Substituting into (0), we obtain an expression for 
tr £} as a sum over all prime cycles p and their repetitions 

trC*- VT V gftzlgg) 
~r P ^|det(l-J;)r 

A Laplace transform replaces the above sum of Dirac 
delta functions by the trace formula for classical 
flows O: 



tr£(s) 
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(10) 

We should caution the reader that in taking the Laplace 
transform we have ignored a possible t — > 0+ volume 
term, as we do not know how to regularize the delta 
function kernel in this limit. In the quantum (or heat 
kernel) case this limit gives rise to the Weyl or Thomas- 
Fermi mean eigenvalue spacing. A more careful treat- 
ment might assign such volume term some interesting 
role in the theory of classical resonance spectra. 

The semi-classical evaluation of the quantum trace is 
consi derab ly more laborious, but the final result, given in 
sect. IV A, is very similar in form to the above classical 
trace. 



C. Fredholm determinants 

The problem with the classical (^), ( [Tof ) and the 
Gutzwiller trace formulas (^0|) is that they diverge pre- 
cisely where one would like to use them (we return to 



this in sect. [VF). While in the physics literature on 
dynamically generated strange sets this does not prevent 
numerical extraction of reasonable "thermodynamic" av- 
erages, in the case of the Gutzwiller trace formula this 
leads to the perplexing observation that crude estimates 
of the radius of convergence seem to put the entire phys- 
ical spectrum out of reach. This problem is cured by 
going from trace formulas to determinants, which turn 
out to have larger analyticity domains. For maps, the 
two are related by 
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F{z) = det(l - zC) = exp - 2^ ~ tr£ " 



For flows the classical Fredholm determinant is given by 
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and the classical trace formula (|io| ) is the logarithmic 
derivative of the classical Fredholm determinant 



(12) 



trC(s) = — logF(s) 



With z set to z — e s , the Fredholm determinant ( |11| ) ap- 
plies both to maps and flows. A Fredholm determinant 
can be rewritten as an infinite product over periodic or- 
bits, by noting that the r sum in jTl| ) is close in form 
to expansion of a logarithm. We cast it into such form 
by expanding the Jacobian weights in terms of stability 
eigenvalues. For a 3-dimensional Hamiltonian flow with 
one expanding eigenvalue A, and one contracting eigen- 
value 1/A, the weight in ([ll]) may be expanded as follows: 

1 1 
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With this we can rewrite the Fredholm determinant ex- 
ponent as 



y> 1 e s7 > 
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and represent the Fredholm determinant as a Selbcrg- 
type product jl7| 
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z is a book-keeping variable that we will use to expand 
zeta functions and determinants, set to z = 1 in calcula- 
tions. In general, t p depends on the dynamical average 
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one wishes to evaluate; this particular weight is used in 
evaluation of escape rates and correlation spectra [fj 18 1. 



The above heuristic manipulations are potentially dan- 
gerous, as we are dealing with infinite-dimensional vector 
spaces and singular integral kernels; we outline now the 
ingredients of the proofs that put the above formulas on 
solid mathematical footing. 



III. FREDHOLM DETERMINANTS CAN BE 
ENTIRE 

As the introduction of the quantum Fredholm determi- 
nant (0) is motivated by its close analogy with the classi- 
cal Fredholm determinants, we shall sketch here the basic 
ideas behind the proofs that the classical Fredholm de- 
terminants are entire, without burdening the reader with 
too many technical details (rigorous treatment is given 
in refs. ||lT| , |l2| , |l3f ) . The reason why one cares whether a 
Fredholm determinant is entire or not is that in practice 
one can extract many more eigenvalues, and to a higher 
accuracy, from entire Fredholm determinants than from 
functions which are not entire. The main point of the 
theorems explained below is that the Fredholm deter- 
minants are entire functions in any dimension, provided 
that 

1. the evolution operator is multiplicative along the flow, 

2. the symbolic dynamics is a finite subshift, 

3. all cycle eigenvalues are hyperbolic (sufficiently bound- 
ed away from 1), 

4. the map (or the flow) is real analytic, ie. it has a 
piecewise analytic continuation to a complex extension 
of the phase space. 

The notion of Axiom A systems is a mathematical ab- 
straction of 2 and 3. It would take us too far to give 
and explain the definition of the Axiom A systems (see 
refs. |]l5|j23|| ) . Axiom A implies, however, the existence of 
a Markov partition of the phase space (see below) from 
which 2 and 3 follow. Properties 1 and 2 enable us to 
represent the evolution operator as a matrix in an ap- 
propriate basis space; properties 3 and 4 enable us to 
bound the size of the matrix elements and control the 
eigenvalues. To see what can go wrong consider the fol- 
lowing examples: 

Property 1 is violated for flows in 3 or more dimensions 
by the following weighted evolution operator 



C t (y,x) = \A t (x)f6(y-f(x)) 



(15) 



where A*(x) is an eigenvalue of the Jacobian transverse 
to the flow. While for the Jacobians J a b — J a Jb for two 
successive segments a and b along the trajectory, the cor- 
responding eigenvalues are in general not multiplicative, 
A a i, ^ A a A& (unless a, b are repeats of the same cycle, so 
JaJb = J ria+nb ) 1 so the above evolution operator is not 



multiplicative along the trajectory. The theorems require 
that the evolution be represented as a matrix in appro- 
priate polynomial basis, and thus cannot be applied to 
non-multiplicative kernels, ie. kernels that do not satisfy 
the semi-group property £} £* = C f +t . 
Property 2 is violated by the 1-d map 

f(x) = a(l - |1 - 2x\) , 1/2 < a < 1 . 

All cycle eigenvalues are hyperbolic, but the critical point 
x c — ^ is in general not a pre-periodic point, there is 
no finite Markov partition, the symbolic dynamics does 
not have a finite grammar, and the theorems discussed 
below do not apply. In practice this means that while 
the leading eigenvalue of C might be computable, the 
reminder of the spectrum is very hard to control; as the 
parameter a is varied, nonleading zeros of the Fredholm 
determinant move wildly about. 
Property 3 is violated by the map 
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Here the interval [0, 1] has a Markov partition into the 
two subintcrvals Iq and I± on which / is monotone. How- 
ever, the fixed point at x — has stability A = 1, and 
violates the condition 3. This type of map is called in- 
termittent and necessitates much extra work [ [l9|| . The 
problem is that the dynamics in the neighborhood of a 
marginal fixed point is very slow, with correlations de- 
caying as power laws rather than exponentially. 

The property 4 is required as from a mathematical 
point of view, the heuristic approach of sect. |h| faces two 
major hurdles: 

1. The trace ([!]) is not well defined since the integral 
kernel is singular. 

2. The existence and properties of eigenvalues are by 
no means clear. 

Both problems are related to how one defines the func- 
tion space on which the evolution operator acts. As 
in physical applications one studies smooth dynamical 
observables, we restrict the space to smooth functions, 
more precisely, the space of functions analytic in a given 
complex domain and having a continuous extension to 
the boundary of the domain. In practice "real analytic" 
means that all expansions are polynomial expansions. In 
order to illustrate how this works in practice, we first 
work out a few simple examples. 



A. Expanding maps 

We start with the trivial example of a repeller with 
only one expanding linear branch 



f(x) = Ax \A\ > 1 
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The action of the associated evolution operator is 

&f>(y) = J dxS(y - Ax)4>(x) = -^-My/A) . 

From this one immediately identifies eigenfunctions and 
eigenvalues: 



1 
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0,1,2,. 



(16) 



The ergodic theory, as presented by Sinai [g0[ and others, 
tempts one to use a space of either integrable or square 
integrable functions. For our purposes, this space is too 
big; had we not insisted on analyticity, non-integer and 
even complex powers could be used in the construction. 
In particular, in the space L 1 all |A|" _1 with a complex 
but Re(a) < 1 would be eigenvalues, ie. the eigenvalues 
would fill out the unit disk. We note that the eigenvalues 
A~" _1 fall off exponentially with n, and that the trace 
of C is given by 



1 OO j 

txC= ]a\ E A ~" = lAKl-A-i 
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in agreement with (^). A similar result is easily obtained 
for powers of £, and for the Fredholm determinant one 
obtains: 



det(l - zC) = Y[(1 



k=0 



|A|A 
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t = —z/\A\, where the coefficients cj, are given explicitly 
by the Euler formula pl| 
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Cfc = 
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(17) 



The coefficients decay asymptotically faster than expo- 
nentially, as A -fc ( fc-1 )/ 2 . This property ensures that for a 
repeller consisting of a single repelling point the classical 
Fredholm determinant is entire in the complex z plane. 

While it is not at all obvious that what is true for a 
single fixed point should also apply to a Cantor set of 
periodic points, the same asymptotic decay of expansion 
coefficients is obtained when several expanding branches 
are involved. Consider a monotone and expanding 1- 
d map f(x), with |/ e '(x)| > 1 on two non-overlapping 
intervals 



,' fo(x), x e h 

~\fi(x), xeh 



(18) 



The simplest non-trivial example is a piecewise-linear 2- 
branch repeller with slopes Ao and Ai. By the chain rule 
A p = Aq^A™ 1 , where the cycle p contains no symbols 
and n\ symbols 1, so the trace (13) reduces to 



tr£" = V i m \ - 



E 



1 



1 



fe - V|A |Ag \A X \A\, 
The Fredholm determinant ( |ll| ) is given by 

where t e = z/|A e |. The eigenvalues (compare with (jl 
are simply 

1 1 



|A |Ag |Ar|Af ' 

Asymptotically the spectrum is dominated by the lesser 
of the two fixed point slopes A — A (if | A 1 < | A x | , other- 
wise A = Ai), and the eigenvalues e~ Vk fall off exponen- 
tially as 1/A fc , just as in the single fixed-point example 
above. 

The proof that the Fredholm determinant for a general 
nonlinear 1-d map ( [l8| ) is entire uses the expansion 

det(l - zC) = ^{-z) k tr (A k C) 

k>0 

where A k C is the k'th exterior power of the operator C. 
For example, 



A 2 £(xia;2,yiy2) = — { 



£(a;i,yi) £(x 2 ,yi) 
C(xi,y 2 ) C{x2,y2) 



so Tr (A 2 C) = i ((Tr Cf - Tr (£ 2 )) . In a suitable poly- 
nomial basis <j) n (z) the operator has an explicit matrix 
representation 



(£</»)„ (z)= T C 



m— 

In the single fixed-point example (fl6|), cf> n = y n , and C is 
diagonal, C nn = A~ n /|A|. 

The proof proceeds by employing Cauchy complex con- 
tour integrals in order to verify that the traces are indeed 
given by the heuristic formula (Q). Furthermore, from 
bounds on the elements C mn one calculates bounds on 
tr (A k £) and verifies |I|@f3| that they again fall off as 

A~ k ~/ 2 , concluding that the C eigenvalues fall off expo- 
nentially for a general Axiom A 1-d map. The simplest 
example of how the Cauchy formula is employed is pro- 
vided by a nonlinear inverse map ip — f , s — sgn(?//) 



C(j>{w) 



dx S(w — f{x))4>{x) = s ip'(w) <j)(ip(w)) . 



Assume that ip is a contraction of the unit disk, i.e. 
\ip{w)\ < 9 < 1 and \ip'(w)\ < C < oo for \w\ < 1 , 
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and expand in a polynomial basis by means of the 
Cauchy formula 



n>0 



dw 4>{w) 



2ni w 



dw <p(i 



2iri w n+1 



In this basis, £ is a represented by the matrix 

dw s i[)'(w)(ip(w)) T 



C(j){w) 



2ni 



„m+l 



to the space of analytic functions, and assume that there 
exists a set of complex neighborhoods Z), D I L such that 
ipij : Cl(Dj) — > Int(-Dj). Mapping closures of domains 
into interiors is a useful way of stating the contraction 
property. The result is pT| , ^2|j2^Jl^] that the expansion 

coefficients Ck fall off as 1/ A h , and the eigenvalues 

as gain, the results can be proved using 

a multinomial basis on each domain and deducing from 
this an explicit matrix representation for the operator. 



Taking the trace and summing we get: 

dw s ip'(w) 



Tr C — ^ ' L nn — 



2ni w — ip(w) 



This integral has but one simple pole at the unique fix 
point w* = ip(w*) = f(w*). Hence 



Tr C 



sip'(w*) _ 1 

i-i>'(w*) ~ \f'(w*) - 1| 



The requirement that map be analytic is needed to guar- 
antee the inequality 

\L mn \ < sup \i/>(w)\ n < C6 n 

\w\<l 

which shows that finite [N x N] matrix truncations ap- 
proximate the operator within an error exponentially 
small in N. 

We note in passing that for 1-d repellers a diagonal- 
ization of an explicit truncated C mn matrix yields many 
more eigenvalues than the cycle expansions The 
reasons why one persists anyway in using the periodic 
orbit theory are partially aestethic, and partially prag- 
matic. Explicit C mn demands explicit choice of a basis 
and is thus non- invariant, in contrast to cycle expansions 
which utilize only the invariant information about the 
flow. In addition, we do not know how to construct C mn 
for a realistic flow, such as the 3-disk problem, while the 
periodic orbit formulas are general and straightforward 
to apply. 

It is a relatively simple task to generalize the above ar- 
guments to an expanding d-dimensional dynamical sys- 
tem / : M — > M with the Markov property (for a more 
precise definition see ref. 

Markov property: One demands that M can be divided 
into S subsets {Iq,Ii, • ■ ■ , Is-x] such that either /Ij n 
Ij = or Ij C fit- The transition matrix takes values 
tij = or 1, accordingly. 

Expansion property: Each inverse ipij : Ij — > U (de- 
fined when — 1) is unique and a contraction. 

Depending on the smoothness of the functions ipij and 
the function space considered, as well as how one defines 
the contraction, one obtains stronger or weaker results on 
the spectrum of eigenvalues. We shall restrict ourselves 



B. Hyperbolic maps 

Theorem [Rugh 1992]: The Fredholm determinant 
for hyperbolic analytic maps is entire. 

The proof, apart from the Markov property which is 
the same as for the purely expanding case, relies heav- 
ily on analyticity of the map in the explicit construction 
of the function space and its basis. The basic idea of 
the proof is to view the hyperbolicity as a cross product 
of a contracting map in the forward time and another 
contracting map in the backward time. In this case the 
Markov property introduced above has to be elaborated 
a bit. Instead of dividing the phase space into intervals, 
one divides it into rectangles. The rectangles should be 
viewed as direct product of intervals (say horizontal and 
vertical), such that the forward map is contracting in, 
for ex. the horizontal direction, while the inverse map 
is contracting in the vertical direction. For Axiom A 
systems the natural coordinate axes are given by the 
stable/unstable manifolds of the map. With the phase 
space divided into S rectangles {Ro, Ri,..., i?s_i}, Ri = 
il 1 x IV with complex extension D\ x D v , the hyperbol- 
icity condition (which at the same time guarantees the 
Markov property) can be formulated as follows: 

Analytic hyperbolic property: Either fRi flint (Rj) = 
0, or for each pair w h i Cl(D^), z v e C\{DV) 
there exist unique analytic functions of Wh,z v : w v = 
w v (w hl z v ) e Int^), z h = z h (w hl z v ) E Int(D^), such 

that f(wh,w v ) = (zh, z v ). Furthermore, if Wh £ i,- and 
z v e IJ, then w v G I V L and z h G ij 1 . (See fig. 0). 

What this means is that it is possible to replace coor- 
dinates Zh,z v at time n by the contracting pair z^,w v , 
where w v is the contracting coordinate at time n + 1 for 
the inverse map. Specifying the closure/interior of the 
sets is a convenient way of defining hyperbolicity. 

A map / satisfying the above condition is called ana- 
lytic hyperbolic and the theorem states that the associ- 
ated Fredholm determinant is entire, and that the trace 
formula (Q) (derived heuristically in ref. |J) is correct. 
We refer the reader to ref. [jl3| for the details of the proof. 
The theorem applies also to hyperbolic analytic maps in 
d dimensions and smooth hyperbolic analytic flows in 
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(d + 1) dimensions, provided that the flow can be re- 
duced to a map by suspension on a Poincare section com- 
plemented by an analytic "ceiling" function p|| which 
accounts for a variation in the section return times. For 
example, if we take as the ceiling function g t (x) = e sT ^ x \ 
where T(x) is the time of the next Poincare section for a 
trajectory staring at x, we reproduce the flow Fredholm 
determinant (|l4|). 

Examples of analytic hyperbolic maps are provided by 
small analytic perturbations of the cat map (where the 
Markov partitioning is non-trivial |24j|), the 3-disk rc- 
peller, and the 2-d baker's map, the last two examples to 
be discussed further below. 

The proofs of discreteness of the classical spectra have 
so far not been extended to the semi-classical Gutzwiller- 
Voros zeta functions, with exception of the spaces of con- 
stant negative curvature l^sfl . The technical problem is 
that the proofs require an evolution operator that is mul- 
tiplicative along the trajectory; composition of the semi- 
classical Green's functions is not of that type, as every 
composition of successive semi-classical paths requires a 



further saddle-point approximation (see sect. IV A) 



trG(E)=g(E) + -Y,T p J2- 



det (1 - 31 



(20) 



Here T p is the p-cycle period, S p its action, m p the Maslov 
index, and 3 P is the transverse Jacobian of the flow. As in 
many applications the wave number k is a more natural 
choice of variable than the energy E, we shall henceforth 
replace E — > k in all semi-classical formulas. 

For 2-d flows the Gutzwiller trace formula is of the 



form (see (13)) 



txG(k)=g(k) 



1 Vt V- 

ih Z. *^|a;|V2(i-i/a;) 



" "k Sp{k)r-\-iirm p r /2 



(21) 



and the corresponding determinant ( |l2| ) is the 2-d 
Gutzwiller- Voros zeta function 



Z qm {k) = cxp - 




D —%S p +iirm p /2 



t p Z 



(22) 



IV. CLASSICAL AND QUANTUM 
DETERMINANTS 

Though the only objective of this paper is to com- 
pare the convergence of the semi-classical determinants 
(|l|) and (|J), theoretical motivation demands a plethora of 
related zeta functions and determinants. In this section 
we introduce and define the additional determinants and 
zeta functions that will be required in what follows. 



A. The Gutzwiller- Voros zeta function 

The semi-classical periodic orbit theory for hyperbolic 
flows was developed by Gutzwiller in terms of traces 
of Van Vleck semi-classical Green's functions ||, and 
subsequently re-expressed in terms of determinants by 
Voros j|. In contrast to the sharp delta- function kernel 
for the classical evolution, the quantum evolution opera- 
tor is the smeared-out Green's function kernel 



G(q",q';t) 



(19) 



defined only on a half of the phase space (typically ei- 
ther the spatial coordinates q, or the momentum coordi- 
nates p). Gutzwiller trace formula follows by replacing 
the quantum Green's function by the semi-classical Van 
Vleck propagator K c (q" ,q';t), and evaluating the trace 
tr K c (q" , q'\ t) by saddle-point methods. The result is a 
periodic orbit formula of topologically the same structure 
as the classical trace (]l^), but with different weights: the 
semi-classical Gutzwiller trace formula B : 



z is a book-keeping variable that keeps track of the topo- 
logical cycle length n p , used to expand zeta functions 
and determinants (see sect. [v]). Unlike the classical Fred- 
holm determinant, which is exact, the Gutzwiller- Voros 
zeta function is the leading term of a semi-classical ap- 
proximation, and the size of corrections to it remain un- 
known |2if |. 

The Gutzwiller trace formula, apart from the quantum 
and Maslov phases, differs from classical trace formula in 
two aspects. One is the volume term g(E) in (^) which 
is missing from our version of the classical trace formula. 
While an overall pre-factor does not affect the location of 
zeros of the determinants, it might play a role in relations 
such as functional equations for zeta functions. The other 
difference is that the quantum kernel leads to a square 
root of the cycle Jacobian determinant, a reflection of the 
relation probability = amplitude 2 . The l/^/dct(l — J p ) 
weight leads in turn to the product representation (0) 



^n(fc) = n n I 1 - wa£) 

p k=0 



(23) 



which differs from the classical Fredholm determinant 
( |i"4| ) by missing exponent k + 1. 



B. The quantum Fredholm determinant 

The conjecture tested in this paper asserts that one 
may replace the Gutzwiller- Voros zeta function (^3|) by 
the quantum Fredholm determinant (|l|) , ie. the Fredholm 
determinant ( |l4| ) with the quantum weights t pi with- 
out disturbing the leading semi-classical eigenvalues, but 
improving the convergence of cycle expansions used in 
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evaluating the spectrum, and revealing part of spectra 
inaccessible to the Gutzwiller-Voros zeta function (see 



sect. VIA for an example). 

The form of the quantum weight (^2|) suggests that the 
quantum evolution operator should be approximated by 
a classical evolution operator with a quantum weight: 



£*(!/, x) = 8(y - f(x))yW&)\e~ iS ^+"™^>/ 2 . 



As explained in sect. [II, this operator is not multiplica- 
tive along the trajectory, and consequently does not sat- 
isfy the assumptions required by the theorems that guar- 
antee that a Fredholm determinant is entire. Neverthe- 
less, our numerical results support the conjecture that the 
Al 1 / 2 weighted determinant has a larger domain of an- 
alyticity than the commonly used Gutzwiller-Voros zeta 
function and that some related determinant might even 
be entire. 



C. The dynamical zeta function 



by 



The Ruelle or the dynamical zeta function M is defined 



1/C = exp 



p r—1 



and has the Euler product representation 



vc=n 



(i - t P ) 



(24) 



(25) 



transport of fc-forms transverse to the flow is given by 
the exterior products of the Jacobian (||) of the flow 

C\ k) {y,x)=8{y-f{x))-A k i. (26) 

In d dimensions (transverse to the flow), the non- 
vanishing exterior powers are J, A X J, A 2 J, . . ., A d J. 
The k = case describes scalar density transport (||); 
k = 1 describes vector transport (utilized in evaluation of 
strange set stabilities ||,|26|,|| and fast dynamo rates p7J ) ; 
k = 2 describes dxAdy area transport, and so on. £^is 

a [d k xd k ] matrix; for example, in the explicit index nota- 
tion A 2 J is given by the antisymmetric exterior product 

(A 2 J) QC , M = \ {8181 ~ 8{8t) J eb J fd . 

As (A fc A) (A fc B) = A k AB, £l n are multiplicative oper- 
ators along the flow, so if det(l — C) is entire, det(l — C^) 
are also entire. The trace formula for Ciia is 



~ tr A fc JI 
tr C (k) (a) = V T„ V , V . pJ 



(fc) 



(27) 



where tr J = Ai+A 2 + . . . + A d , tr A 2 J = A1A2 + A1A3 + 
. . . + A^-iArf, and so on. The corresponding Fredholm 
determinants are: 
1-d case: 



det(l-£ (1) )=niI( 1 -WAp 

p fc=0 



k-l\ 



where the product is over all prime cycles p, and the cycle 
weight t p depends on the average computed; the simplest 
example is the weight (|l4| ) used in computation of escape 
rates and correlation spectra. ( |25| ) also yields the leading 
semi-classical quantum resonances, if t p is the quantum 
weight (|2^) associated with the cycle p. 

Historical antecedents of the dynamical zeta func- 
tion are the fixed-point counting functions introduced 
by Artin-Mazur and Smale p5|, and the determinants 
of transfer operators of statistical mechanics |20|. While 



l/£ is the natural object in these applications, in dy- 
namical systems theory dynamical zeta functions arise 
naturally only for piecewise linear mappings; for smooth 
flows the natural object for study of classical and quan- 
tal spectra are the determinants introduced above. How- 
ever, dynamical zeta functions will here be useful objects 
in relating various determinants and explaining the geo- 
metrical meaning of curvature expansions. 



D. Weighted Fredholm determinants 



The theorems of sect. [II apply not only to the evo- 
lution operator (||), but also to more general evolution 
operators multiplicative along the flow. In particular, 



2-d Hamiltonian case 
det(l — £(i) 



nn<i 

p k=0 



A 



k-1 



dct(i-£ (2) )=nn( i -wA^) 

p k=0 



k+1 



k+l 



A k+1 



k+l 



det(l - C) , 



where we have used the Hamiltonian volume conservation 
in tr A 2 J p = A p • 1/A p = 1. 

In order to discuss relations between various determi- 
nants and zeta functions, it is convenient to also define 
following weighted dynamical zeta functions: 



1 tl 



r A kr 

p r=l P 



(28) 



and weighted Fredholm determinants 

i (WApT 



F, 



r \2 



(29) 



(Fq = F) . For 2-dimensional Hamiltonian flows the clas- 
sical Fredholm determinant (n4h can be written as 
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F = 



nvc 



fc+i 

k 



(30) 



fc=0 



Fk can be interpreted as the Fredholm determinant 
det(l — Ck) of the weighted evolution operator 
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A t (x)' 



where A' (x) is the expanding eigenvalue of the Jacobian 
transverse to the flow, and g t {x) is any analytic weight 
function multiplicative along the trajectory. As discussed 
in sect. [II, there is no guarantee that det(l — Ck), the 
Fredholm determinant for this operator, should be entire 
for k ^ 0. 



E. Relations between determinants and zeta 
functions 

Ruellc has observed that the elementary identity 
for d-dimensional matrices 



1 



det(l - J) 



£(-l) fe tr (A fe J) , 



fe=0 



inserted into the exponential representation (|24) of the 
dynamical zeta function, relates the dynamical zeta func- 
tion to weighted Fredholm determinants. In one dimen- 
sion this identity 



1 



1 



1 



1 



(1 - 1/A) A (1 - 1/A) 



expresses the dynamical zeta function as a ratio of two 
Fredholm determinants 



i/C 



det(l - £) 
det(l-£ (1) ) 



(31) 



and shows that 1/C is meromorphic, with poles given by 



the zeros of det(l — £(i))- In refs. [p|,|26| heuristic ar- 
guments were developed for 1-dimensional mappings to 
explain how the poles of individual 1/Cfe cancel against 
the zeros of 1/Cfe+i, and thus conspire to make the cor- 
responding Fredholm determinants entire. Numerical 
checks verify both the heuristic arguments, and the for- 
mula di]). 

For 2-dimensional Hamiltonian flows this yields 



i/C 



det(l-£)det(l -£ (2) ) _ F 2 
det(l-£ (1) ) = F\F~i 



This establishes that 1/C, is meromorphic in 2-d as well, 
but the relation is not particularly useful for our pur- 
poses. Instead we insert the identity 



1 = 



1 



1 



(1-1/A) 2 A(l-1/A) 2 A 2 (l-1/A) 2 



into the exponential representation 
tain 



i/C* 



FkFk 



+ 2 



F 2 

r k+l 



s) of 1/Oc, and ob- 



(32) 



Even though we have no guarantee that Fk are entire, 
we do know (by arguments explained in the next sec- 
tion) that the upper bound on the leading zeros of -Ffc+i 
lies strictly below the leading zeros of Fk, and therefore 
we expect that for 2-dimensional Hamiltonian flows the 
dynamical zeta function l/0fe has generically a double 
leading pole coinciding with the leading zero of the Fk+i 
Fredholm determinant. This might fail if the poles and 
leading eigenvalues come in wrong order, but we have 
not encountered such situation in our numerical investi- 
gations. This result can also be stated as follows: the 
theorem that establishes that the classical Fredholm de- 
terminant d3(j| ) is entire, implies that the poles in 1/Oc 
must have right multiplicities in order that they be can- 
celled in the F = J] 1/Cfc +1 product. 

Both the quantum Fredholm determinant and the 
quantum zeta function yield the same leading zeros, given 
by 1/Co- They differ in non-leading zeros (zeros with 
larger imaginary part of the complex energy), but as the 
Gutzwiller-Voros zeta function (g3|) is only the leading 
term of a semi-classical approximation, with the size of 
corrections unknown, the physical significance of these 
non-leading zeros remains unclear. 



F. Abscissa of absolute convergence 

Consider the "thermodynamic" approximation 
the case of the Gutzwiller trace formula (Ell): 



trG(fc) « E T pE' 

p r 



r S p ( k ) r + i 7T m p r / 2 



\K\ 1/2 



This approximation omits the non-leading « 1 / A p terms 
that vanish in the t — > oo limit and do not affect the 
leading eigenvalues. If the phases conspire to partially 
cancel contributing terms, the sum diverges for a larger 
value of Im(fc). The abscissa of absolute convergence in 
the complex k plane is obtained by maximizing the sum, 
ie. replacing all terms by their absolute values: 



s T p Im(fc)r 



(we have for simplicity taken S p to be the action for bil- 
liards, Sp/h = T p k). Value of Im(fc) for which this sum 
diverges determines the abscissa of absolute convergence. 
In the case of the Gutzwiller trace formula (|2l| ) this leads 
to the disappointing result that all eigenvalues of interest 
are outside the domain of convergence. However, for de- 
terminants and zeta functions the eigenvalues are given 







by locations of the zeros, and the analyticity domain is 
larger. We can also use determinants to accurately esti- 
mate this abscissa of absolute convergence, by replacing 
all cycle weights t p in (O) by their absolute values. 

To evaluate the abscissa of absolute convergence of the 
Gutzwiller-Voros zeta function we first note that insert- 
ing the identity 1 = (1 - 1 /A£)/(l - 1 /Ap into the expo- 
nent of Gutzwiller-Voros zeta function ( |22| ) , one obtains 
the following relation between the Gutzwiller-Voros zeta 
function and the quantum Fredholm determinant: 



F qm {k) 



where 



Fi (fc) = exp 



rA;(l-l/A;) 2 



(33) 



The radius of convergence of Z qm (k) is therefore deter- 
mined by the leading zeros of Fi(k). To estimate the 
upper bound on Im(&;) for the zeros of Fi , we omit all 
signs and phases in the weights in Fi : 



nW = -P^-EE;| A;Kl _ 1/A;)S 



(34) 



and compute it s lead ing zero at Re(fc) = 0. An example 
is given in sect. VIA. 



G. The Gutzwiller-Voros zeta function puzzle: a 
classical determinant? 

On the basis of close analogy between the classical and 
the quantum zeta functions, it has been hoped J^] that 
for Axiom A systems the Gutzwiller-Voros zeta functions 
( ^3| ) should also be entire. This hope was dashed by Eck- 
hardt and Russberg |3l| ] who have established by numer- 
ical studies that the Gutzwiller-Voros zeta functions for 
3-disk repeller have poles. They find numerically that 
in the Gutzwiller-Voros zeta function (|23|) product rep- 
resentation 1/Co has a double pole coinciding with the 
leading zero of 1/Ci- Consequently l/£o, 1/CoCi an d Z qm 
all have the same leading pole, and coefficients in their cy- 
cle expansions fall off exponentially with the same slope . 
Our numerical tests on the 3-disk system (see sect. VIA ) 
support this conclusion. 

Why should l/£o have a double leading pole? The dou- 
ble pole is not as surprising as it might seem at the first 
glance; indeed, for the classical Fredholm determinant 
( |l4| ) 1/Co must have a double pole in order to cancel 
the leading double zero of the (1/Ci) 2 factor. In other 
words, the numerics indicates that the semi-classical de- 
terminants and zeta functions have analyticity proper- 
ties characteristic of classical determinants. As a priori 
any determinant or zeta function which yields the same 



leading semi-classical resonances is equally good, one is 
lead to pose the conjecture tested in this paper: quan- 
tum Fredholm determinants are preferable to Gutzwiller- 
Voros zeta functions, as they yield the same leading spec- 
trum, but with better analyticity properties. 

As we have no proof, we resort to numerical tests of 
the conjecture; but first we need to briefly explain how 
such determinants are investigated numerically. 



V. CYCLE EXPANSIONS 

A cycle expansion [flo| is a series representation of a 
zeta function or a Fredholm determinant, expanded as a 
sum over pseudo-cycles, products of prime cycle weights 
t p , ordered by increasing cycle length and instability. The 
products ((2^), ((2^) are really only a shorthand notation 
for zeta functions and determinants - for example, the 
zeros of the individual factors in infinite products j2^), 
( p5| ) are not the zeros of the corresponding zeta functions 
and determinants, and convergence of such objects is far 
from obvious. 



A. Curvature expansions 

Curvature expansions are based on the observation 
ifjojpll that the motion in dynamical systems with finite 
grammar is organized around a few fundamental cycles; 
more precisely, that the cycle expansion of the dynami- 
cal zeta function ( p5| ) allows a regrouping of terms into 
dominant fundamental contributions tf and decreasing 
curvature corrections c„: 



i/c = i-£*/-£. 



(35) 



/ 



The fundamental cycles tf have no shorter approximants; 
they are the "building blocks" of the dynamics in the 
sense that all longer orbits can be approximately pieced 
together from them. In piecewise linear approximations 
to the flow, 1/C is given by the determinant for a finite 
Markov transition matrix, and all c„ vanish identically. 
Hence the designation "curvatures" ; size of c n is an indi- 
cation how far the flow is from a piecewise linearization. 

A typical curvature term in ([55]) is a difference of a long 
cycle {ab} and its shadowing approximation by shorter 
cycles {a} and {&}: 

tab — t a tb — t a b(l — t a tb/t a b) 

The orbits that follow the same symbolic dynamics, such 
as {ab} and the "pseudo orbit" {a} {&}, lie close to each 
other, have similar weights, and for longer and longer 
orbits the differences are expected to fall off rapidly. 
For systems that satisfy Axiom A requirements, such as 
the 3-disk repeller, curvature expansions converge very 
well Jig]. It is crucial that the curvature expansion is 
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grouped (and truncated) by topologically related cycles 
and pseudo-cycles; truncations that ignore topology, such 
as inclusion of all cycles with T p < T max , will contain 
un-shadowed orbits, and exhibit a mediocre convergence 
compared with the curvature expansions. 



Russberg and Eckhardt |31|] have generalized curva- 
ture expansions to Fredholm determinants; they work to 
a point, but are less powerful than the general theorems 
on Axiom A flows. For example, they predict that the 
classical Fredholm determinant could have poles, in vio- 
lation of the theorem explained above. 



B. Fredholm determinant cycle expansions 

While for the dynamical zeta function cycle expansions 
the shadowing is easy to explain, the resulting conver- 
gence is not the best achievable; as explained above, 
Fredholm determinants are expected to be entire, and 
their cycle expansions should converge faster than expo- 
nentially. The Fredholm determinant cycle expansions 
are somewhat more complicated than those for the dy- 
namical zeta functions. We expand the exponential rep- 
resentation ( pi] ) of F(s) as a multinomial in prime cycle 
weights t p 



oo 

= E c P kt p 



v 1 

, =1 r|det [l-J;]| + 2 ( - 



This yields the cycle expansion for F(s): 



F(s) 



T k-t fco kr> 

Pi p 2 p 3 ■ 



E 

oo 



T k-t feo kr> 
Pi >2 ?>3 



where the sum goes over all pseudo-cycles. The coeffi- 
cients have a simple form only in l-d, given by the Euler 
formula (|l7|). Expansions for the 2-d case are discussed 
in refs. J22| , pl[ . 

In practice we do not do anything as complicated. We 
evaluate numerically the exponent in ( p"l| ) as a power se- 
ries in the book-keeping variable z, truncated to maximal 
cycle length N. The we expand the exponential, keeping 
terms up to N and obtain 



F N (s,z) 



N 

E 

fe=0 



C k (s)z k 



(36) 



C. Convergence of cycle expansions 

It is fairly easy to establish that for Axiom A systems 
the trace formulas converge exponentially with the num- 
ber of cycles included. As explained in sect. [VF , the 
trace formulas are not absolutely convergent where you 
need them, and in addition, shadowing of longer orbits 
by nearby pseudo-orbits is not implemented, so we will 
not use trace formulas at all. However, it should be noted 
that for systems other than Axiom A, we do not know 
how to improve convergence by shadowing cancellations, 
or define determinants that are guaranteed to be entire, 
and it is still possible that for generic systems determi- 
nants do not converge any better than traces. 

For dynamical zeta functions geometrical estimates || 
imply that for Axiom A systems the curvature expan- 
sion coefficients fall off exponentially, C k ~ C k , and the 
expansion sums up to a pole 



E'" 

n=0 



oo 



C k z K ~^{Cz) 

n=0 



1 



1 - Cz 



Such poles are expected from Ruelle's relation (^Tj) be- 
tween dynamical zeta functions and Fredholm determi- 
nants, and are indeed observed numerically ]26] ]. Conver- 
gence of dynamical zeta functions cycle expansions can 
be accelerated by a variety of numerical methods, but 
both on theoretical grounds and in practice, the preferred 
alternative is to use Fredholm determinants instead. 
As shown in sect. 



[II 



for Axiom A maps the coef- 
ficients in the classical Fredholm determinant F(z) = 
^nCnZ" expansions fall off faster than exponentially, 

as C„ ps A~" 2 for l-d maps, as C„ ps A~"' V2 for 2- 
dimensional maps / 3-dimensional flows, and in d dimen- 
sions as A - ™ + . These estimates are confirmed by the 
numerical tests of ref. Jl3| , numerical results for the 3- 
disk repeller ref. J31J1 , as well as the numerical results 
presented below. 



and similarly for the Gutzwiller-Voros zeta function and 
the dynamical zeta functions. In the final evaluation z 
is set to z = 1, but the organization by powers of z k 
is crucial to the convergence of cycle expansions. For 
example, for l-d Axiom A mappings C k (s) « C~ k for 
any s - this super-exponential convergence is precisely 
the reason why the variable z was introduced in the first 
place. The zeros of Fn(s, 1) are determined by standard 
methods, such as the Newton-Raphson algorithm. 



VI. NUMERICAL RESULTS 

In this section we present the evidence that the quan- 
tum Fredholm determinant is numerically as convergent 
as the classical Fredholm determinant, in contrast to the 
Gutzwiller-Voros zeta function which has a finite radius 
of convergence. 
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A. 3-disk resonances 



Following refs. |pL8| , P l|j2S|j30f1 , we start by performing 
our numerical tests on the 3-disk repeller. The 3-disk 
repeller is the simplest physical realization of an Axiom 
A system, particularly convenient for numerical investi- 
gations. The methods that we use to extract periodic 
orbits, their periods and their stability eigenvalues are 
described in ref. For billiards the cycle weight t p 

required for evaluation of the classical escape rates and 
correlation spectra is given by ( |l4| ) . The action S p is pro- 
portional to the cycle period T p , and the Maslov index 
changes by +2 for each disk bounce, m v — 2n p , so the 
quantum weight ((2^) is given by 



t P = (-1)' 



(37) 



where k = (momcntum)/27r is the wave-number, and we 
take velocity = 1, mass = 1. 

Cycle expansion ([56]) coefficients \C n \ for different de- 
terminants and zeta functions are plotted in figs. ^ and |3| 
as function of the topological cycle length n. Zeta func- 
tions exhibit exponential falloff, implying a pole in the 
complex plane, while both the classical and the quantum 
Fredholm determinants appear to exhibit a faster than 
exponential falloff, with no indication of a finite radius 
of convergence within the numerical validity of our cycle 
expansion truncations. 

In particular, the quantum Fredholm determinant en- 
ables us to uncover a larger part of the quantum reso- 
nances than what was hitherto accessible by means of 
the dynamical zeta functions [pO 31 32|. The eye is con- 



veniently guided to the zeros by means of complex s 
plane contour plots, with different intervals of the ab- 
solute value of the function under investigation assigned 
different colors; zeros emerge as centers of elliptic neigh- 
borhoods of rapidly changing colors. Detailed scans of 
the whole area of the complex s plane under investigation 
and searches for the zeros of classical and quantum Fred- 
holm determinants, fig. ^, reveal complicated patterns 
of resonances, with the classical and the semi-classical 
resonance patterns surprisingly similar. It is known fl§| l 
that the leading semi-classical resonances are very accu- 
rate approximations to the exact quantum resonances; 
the semi-classical resonances further down in the s com- 
plex plane in fig. || have not yet been compared with 
the exact quantum values. It would be of interest to 
check whether also all quantum Fredholm determinant 
resonances correspond to the true quantum resonances. 

An interpretation of the resonance spectrum of the 
classical Fredholm determinant is given in ref. |3C| l , where 
the resonances are related to the oscillations in N(t), the 
number of particles that have not escaped by the time t, 
with the basic frequency ZJ — 2n/T given by the inverse 
of T, the mean flight time between the disks. ZJ yields 
the mean spacing of the resonances along the imaginary 
k axis. In the 3-disk system there are two fundamental 



frequencies, u>o and ui\, determined by the inverse periods 
of the two fundamental cycles and 1. Corresponding 
beat frequency J33] / = 2ir/ (T± — T ) is clearly visible in 
fig. ||. A rough measurement of the period of the beats 
in fig. U yields some 23.7 units along the real k axis, to 
be compared to 2n/(Ti - T ) = 23.4. 

Effect of sub-dominant resonances on the measurable 
spectra are exponentially small, and presumably of lit- 
tle physical interest. We investigate them in detail here 
mostly in order to demonstrate that our determinants 
indeed exhibit better convergence than the Gutzwillcr- 
Voros zeta function. 

Contour plots arc also helpful in comparing the domain 
of convergence of the Fredholm determinant to that of 
the Gutzwiller-Voros zeta function. As can be seen from 
fig. ^, the quantum Fredholm determinant can be contin- 
ued considerably farther down in the complex k plane, in 
contrast to the dynamical zeta function scans such as 
those given in ref. J30|. While the zeta functions clearly 
exhibit a finite radius of convergence, in agreement with 
the arguments of sect. [V, both the classical Fredholm 



determinant and the quantum Fredholm determinant be- 
have as entire functions. We compute the abscissa of 
absolute convergence for the Gutzwiller-Voros zeta func- 
tion by means of (|34]); for the case at hand we obtain 
the leading zero at k ac = 0.0 - % 0.699110157151 . . .. In- 
deed, comparison of the contour plots of fig. ^ shows that 
no feature of the Gutzwiller-Voros zeta function coun- 
ter plot below k c is significant. Interestingly enough, 
the apparent border of Gutzwiller-Voros zeta function 
convergence in fig. || seems to coincide with Re(s) = 0, 
Im(s) = —1.09653395 . . ., the zero obtained from Fi/ 2 (k) 
by removing quantum phases, t p — ► \t p \, but keeping the 
eigenvalue A p sign, in eq. (|3^). 

As discussed above, for l-d systems the pole of C^ 1 
coincides with the leading zero of Ci 1 , an d the result- 
ing product remains finite and has a zero close ||||] 
to the leading (^-f 1 zero. In simple examples, such as 
the symmetric 1-dimensional tent map repeller, the non- 
leading / leading zeros of the classical Fq / F\ maps are 
identical. This suggests that some of the non-leading 
zeros of F are shadows of the C-f 1 zeros and hence 
likely to lie close to the leading zeros of F\. For ex- 
ample, the non-leading resonance k = 0.9915231008 + 
i 12.5163342128 of the classical Fredholm determinant 
(see fig. ^), while distinct from the leading resonance of 
Fi at k = 0.99197582677 + i 12.5029206443798, belongs 
to a family of resonances that all lie very close to the 
leading Fi resonances. We are confident that such res- 
onances are distinct, as their cycle expansions converge 
super-exponentially to all digits listed above, in agree- 
ment with the general theory. 

The F spectrum that is not echoed by the F\ spec- 
trum can be isolated by "subtracting" the F\ spectrum, 
ie. removing those resonances from the F spectrum that 
lie closer than some e to a F\ resonance. The resulting 
spectrum is depicted in fig. ^. 
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What is particularly interesting about this spectrum is 
that with the resonances associated with the Fi spectrum 
removed, families that connect non-leading resonances at 
Re(fc) = with the leading part of the resonance spec- 
trum for larger Re(fc) are clearly visible in fig. |^. This 
makes it rather clear that the evaluation of leading res- 
onances for large Re(fc) requires inclusion of longer cy- 
cles, the same that are required to control the spectrum 
for large negative Im(fc) at Re(fc) = 0. This is the co- 
nundrum of semi-classical cycle expansions; while semi- 
classical intuition implies that the Gutzwiller-Voros zeta 
function should be applicable for large Re(fc), in prac- 
tice the semi-classical cycle expansions work best for the 
bottom of the spectrum. 



B. Henon mapping as a normal form for a flow 

Our next example is a model of a 3-dimensional flow 
x = F(x), x = (xi, X2,Xa). We assume that the flow of 
interest is recurrent, and that given a convenient Poincare 
section coordinatized by coordinate pair (x, y), the flow 
can be described by a 2-dimensional Poincare map 



P : 



x' = h(x,y) 
y' = f2(x,y) 



(38) 



together with the "ceiling" (23| function T(x, y) 
which gives the time of flight to the next section 
for a trajectory starting at {x,y). A trajectory 
(xx,yx), (x 2 ,y2), (2:3,2/3), • ■ ■ induces a sequence of flight 
times Ti, T2, T3, . . ., = T(xk,yk)- This sequence can 
be used to construct an embedding space, with the pair 
(Tfc,Tfc+i) serving as the embedding coordinates. The 
embedding theories |33| imply that the sequence (xk, yk) 
generated by the original map is dynamically equiva- 
lent with the sequence (Tk,Tk+i) (provided that the 
coordinate transformation is non-singular, for example 
T =const is excluded), and that the mapping 



M : (Tfc.Tfc+i) -» (Tfc+i.Tfc+a) 



(39) 



gives us the full invariant information about the dynam- 
ical system. In particular, both the return times T/~ and 
the stability eigenvalues of periodic orbits computed from 
( |39| ) are equal to those of the original map (|38|). While 
the functional form of ([39]) can be complicated in actual 
implementation, we shall use here these ideas only as a 
motivation for the following rather simple model: 

Assume that the Poincare map of the flow can be mod- 
elled by a single analytic map for the whole T range of 
interest: 



Tk.+i = G(Tk,Tk-i). 



(40) 



In general G can be a complicated function, but the es- 
sential properties of a continuous flow can be modelled 
by the Henon map, which we take as a local normal form 



(up to quadratic terms) of the time mapping. In the 
Hamiltonian case, the form of the Henon map is 



T k =T + ax k 
Xk+i = 1 - ax\ - x k - 



(41) 



T is essentially the mean time of flight between Poincare 
sections which we shall set to T = 1, and the parameter 
a allows us to choose a narrow or a broad return time 
distribution around the mean return time. For a — 
the flow reduces to the usual Henon map, with constant 
time between the iterations. For example, the R : a = 
6 : 1 3-disk repeller has a rather thin repeller, and can 
be roughly fit with a « 20 in (|lj) . 

We relate this hypothetic flow to a semi-classical sys- 
tem by interpreting the flight times ([D]) as the lengths 
of segments of billiard trajectories, with the action of the 
periodic orbit p given by 



■rS p (k) = k Tj = k(n p + x P ,i) ■ 



(42) 



i=l 



where x Pt i is the coordinate of the ?-th cycle point of 
prime cycle p, and k is the wave number. For billiards, 
the Maslov index m p is taken equal to for cycles with 
positive stability eigenvalue A p , and 1 for A p negative. 
We take unit velocity, so T p — L p . 

For the complete repeller case (all binary sequences are 
realized), the cycles are evaluated as follows. According 
to (41), the coordinates of a periodic orbit of length n p 
satisfy the equation 



X p ,i-\-l ~\~ X p .i — l 



1 



1, 



(43) 



X Pf n p • In 



with the periodic boundary condition x Pj o 
the complete repeller case, the Henon map is a realiza- 
tion of the Smale horseshoe, and the symbolic dynam- 
ics has a very simple description in terms of the binary 
alphabet e = 0,1, e p ^ = (1 + S p ^)/2, where S Pt i are 
the signs of the corresponding cycle point coordinates, 
S Pt i = siga(x Pt i). We start with a preassigned sign se- 
quence S Pl i, <Sp,2, ■ ■ • , S Pt n p , and a good initial guess for 
t he coordinates x' p i . Using the inverse of the equation 



x 



P>»— 1 



a 



i = l,...,n p (44) 



we converge iteratively, at exponential rate, to the desired 
cycle points x p ^. Given the cycle points, the cycle sta- 
bilities and periods are easily computed. The times and 
the stabilities of the short periodic orbits for the Henon 
repeller at a = 6 are listed in table Q; in actual cal- 
culations we use all prime cycles up to topological length 
n = 12 (when needed, all cycles up to length n = 20 
have been computed) . Once we have constructed a table 
of lengths and stabilities, extraction of the eigenvalues 
proceeds as in the three-disk example discussed above. 
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We consider first the discrete time approximation, with 
the parameter a set equal to zero. Then the lengths of 



the periodic orbits are simply T p 



the spectrum 



in the complex wave-number k plane is periodic along 
the real direction, and it is sufficient to consider the 
< Re(fc) < 2vr strip. Fig. § shows well separated eigen- 
values for both the classical Fredholm determinant and 
the quantum Fredholm determinant, with no hint of a 
border of analyticity. In contrast, a corresponding con- 
tour plot of Gutzwiller-Voros zeta function (not included 
here) shows very clearly the finite radius of convergence, 
similar to what is observed in figs, ^(a) and |l(| 

When we switch from a map to a model of a flow by 
introducing a time spread, a ^ in (|42| ), the spectrum 
is no longer periodic in Re(/c). In fig. [8pwe plot the clas- 
sical Fredholm determinant and the quantum Fredholm 
determinant resonances on the same scale as in fig. 0, 
and in fig. || we plot a wider Re(fc) range. Both classical 
Fredholm determinant and the quantum Fredholm de- 
terminant converge as well as in the a = discrete time 
case (where the theorem of sect. Ill guarantees that at 
least the classical Fredholm determinant is entire), and 
in the Re(fc) direction a quasi-periodic pattern of eigen- 
values can be observed, similar to the 3-disk resonance 
spectrum of fig. |[ 

The unexpected feature of this spectrum is that the 
leading quantum resonance in the < Rc(fc) < 2tt do- 
main has larger imaginary part then the leading reso- 
nances in the sectors 2im < Re(fc) < 2ir(n + 1) n = 
1,2, .... This is a consequence of the partial cancellation 
of the contributions from the two fundamental cycles 
and 1 in the cycle expansion (|22]), due to relative minus 
sign arising from the Maslov phases. A partial cancela- 
tion of the leading terms in the cycle expansion dimin- 
ishes the leading coefficients in the I /Co cycle expansion, 
hence the leading zero has a larger imaginary part. In the 
3-disk example these contributions were positive in the 
A\ subspace, but in the A 2 subspace they partially can- 
cel in the same way as in our example. For the classical 
Fredholm determinant, the weight ([l4]) depends only on 
the absolute magnitude of the stability, so the imaginary 
part of the leading zero of the classical Fredholm determi- 
nant in the domain < Re(fc) < 2-k is smaller then the 
imaginary part of the Gutzwiller-Voros zeta function's 
leading zero. This goes contrary to the usual expecta- 
tion that the quantum escape rate should be slower than 
the classical one p9| . 

This model of a flow could also be applied to smooth 
potentials, with other choices for the Maslov indices, 
and the other forms of the action. We have computed 
the spectra with Maslov indices appropriate to flows in 
smooth potentials, and found qualitatively similar ana- 
lytic properties of the quantum Fredholm determinant. 



C. Baker's map 

As the third and the last illustration of the improved 
convergence obtained by use of the quantum Fredholm 
determinants, we evaluate both the quantum Fredholm 
determinant and Gutzwiller-Voros zeta function for a 
model 2-d map that takes the unit square into itself by 
the following piecewise analytic transformations: 

for < x 2 < 2/5: 

2 2 
fi{xi,x 2 ) = -xi + x 2 (- - x 2 )x 1 (l - Xi) 

O 

5 2 
f 2 (xi,x 2 ) = -x 2 + 2x 2 (- - at 2 )a;i(l - xi) 
2 5 



for 2/5 < x 2 < 1: 

3 7 
fi{xi,x 2 ) = — x x + — 

f 2 (xi,x 2 ) = ^{x 2 - jj) . 



(45) 



The model has a binary symbolic dynamics and no par- 
ticular physical motivation; it was introduced in ref. 
to illustrate numerically the theorem that the classical 
Fredholm determinant for an Axiom A hyperbolic flow 
is entire. We study it here merely as another conve- 
nient model for numerical investigations of the conver- 
gence of the contracting map equivalents of the quan- 
tum Fredholm determinant and the Gutzwiller-Voros 
zeta function. For this discrete time example, we as- 
sume no Maslov indices, and chose the action to be 
equal to the (integer) time. The "Gutzwiller trace for- 
mula" for this map is defined by assigning cycles weight 
det(l — J p )|, with the corresponding "Gutzwiller- 
Voros zeta function" . 

The convergence of the cycle expansions (|3^) for the 
two determinants is illustrated by fig. Starting with 
cycle length n — 13, and up to n = 24, the highest 
cycle length computed, the quantum Fredholm deter- 
minant performs better than the Gutzwiller-Voros zeta 
function. The convergence is also illustrated by the ta- 
ble [n| where the zeros of the two determinants are shown. 
Even though the map is not area preserving and the 
"Gutzwiller-Voros zeta function" is quite artificial, the 
convergence is much the same as for the 3-disk repeller. 

The numerical evidence from baker's map on whether 
the quantum Fredholm determinant and weighted Fred- 
holm determinants are entire is inconclusive; zeros of 
F-i, F qm , i*b, F\ cycle expansion truncations up to 
n-max = 22 have qualitatively the same distribution in 
the complex plane. l/C(/c) and Z qm (k) exhibit somewhat 
worse convergence than the Fredholm determinants. In 
all cases, only a few leading zeros are true zeros; for ex- 
ample, table [lj lists 7 true zeros, the remaining 17 form a 
"wedge" in the complex plane (see fig. O), which slowly 
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drifts toward higher Im(fc) with increased cycle trunca- 
tion length. 

We conclude our discussion of numerics with a word 
of caution: as we do not know how quickly the asymp- 
totics should set in, our numerical results can easily be 
misleading. For example, for a larger disk-disk spacing, 
pre-asymptotic oscillations are visible in fig. 2, and one 
might mistakenly conclude |3l]] from such data that the 
classical Fredholm determinant has a pole. There is no 
substitute for theorems that established that appropri- 
ate determinants are entire; such oscillations make nu- 
merical convergence uncertain already in the simplest 1- 
dimensional repellers |3q]. 



VII. CONCLUSIONS AND SUMMARY 

In conclusion, we have tested numerically the con- 
jecture that the new approximation to the quantum 
determinant, the quantum Fredholm determinant, has 
better analyticity properties than the commonly used 
Gutzwiller-Voros zeta function. Existence of such deter- 
minant suggests a starting approximation to the quan- 
tum propagator different from the usual Van Vleck semi- 
classical propagator. As we do not know an operator 
whose determinant is the quantum Fredholm determi- 
nant, we have no proof that such determinant is entire, 
only numerical evidence that its convergence is superior 
to the Gutzwiller-Voros zeta function. The new deter- 
minant could be of practical utility, as for nice hyper- 
bolic systems its convergence is superior to that of the 
Gutzwiller-Voros zeta functions. 

We were guided here by the Axiom A intuition de- 
veloped by Smale, Ruelle, and others: if the dynamical 
evolution can be cast in terms of an evolution operator 
multiplicative along the flow, if the corresponding map- 
ping (for ex., return map for a Poincare section of the 
flow) is analytic hyperbolic, and if the topology of the 
repeller is given by a finite Markov partition, then the 
Fredholm determinant ( |T^) is entire. An alternative ap- 
proach, inspired by the theory of the Ricmann zeta func- 
tion, is due to M.V. Berry and J. Keating ]35|] . The idea 
is to improve the periodic orbit expansions by impos- 
ing unitarity as a functional equation ansatz. The cycle 
expansions used are the same as the original ones |&^|, 
but the philosophy is quite different; the claim is that the 
optimal estimate for low eigenvalues of classically chaotic 
quantum systems is obtained by taking the real part of 
the cycle expansion of the semi-classical zeta function, 
cut at the appropriate cycle period. 

The real life challenge are generic dynamical flows, 
which fit neither schematization. Unfortunately we know 
of no smooth potential which is both Axiom A, and has 
bound states. Most systems of interest are not of the 
"Axiom A" category; they are neither purely hyperbolic 
nor do they have a simple symbolic dynamics grammar. 



The crucial ingredient for nice analyticity properties of 
zeta functions is existence of finite grammar (coupled 
with uniform hyperbolicity) . From hyperbolic dynam- 
ics point of view, the Riemann zeta is perhaps the worst 
possible example; understanding the symbolic dynamics 
would amount to being able to give a finite grammar 
definition of all primes. Hyperbolic dynamics suggests 
that a generic "chaotic" dynamical system should be ap- 
proached by a sequence of finite grammar approximations 
H , pretty much as a "generic" number is approached by 
a sequence of continued fractions. 

The dynamical systems that we are really interested 
in - for example, smooth bound Hamiltonian potentials - 
are presumably never really chaotic, and it is still unclear 
what intuition is more rewarding: are quantum spectra 
of chaotic dynamics in smooth bound Hamiltonian po- 
tentials more like zeros of Riemann zetas or zeros of dy- 
namical zeta functions? 

Note added in proof: Since completion of this work, 
Vattay et al [[3?J have succeeded in constructing a mul- 
tiplicative evolution operator for semi-classical quantum 
mechanics, with the corresponding determinant 



exp 



E 



1|A P | 



-r/3 p i(S p (k)/H+u p 7r/2)r 



r (1 



i/a;) 2 (i 



1/Af) 



(46) 



This determinant is entire for the Axiom A flows. The 
results presented in this paper remain valid, but the new 
formulation makes it possible to reformulate our conjec- 
tures as theorems. In particular, the quantum Fredholm 
determinant can now be written as a ratio of two Vattay 
determinants 



Fqm (^) 



and as explained in sect. [VF, the abcissa of convergence 



for such ratio is given by the upper bound on the leading 
zero of -F(|, fe). For example, for the 3-disk R : a = 6 : 1 
system we find 



Im(fc QC 



-1.276625955.. 



so the quantum Fredholm determinant can be contin- 
ued almost a factor 2 down in the complex k plane be- 
yond the region of applicability of the Gutzwiller-Voros 
zeta function. This explains why for all practical pur- 
poses the quantum Fredholm determinants behave as en- 
tire functions: only careful numerics reveals the differ- 
ence between the quantum Fredholm determinant and 
the Vattay determinant in higher order terms in cycle 
expansions. The difference is illustrated by fig. |1| the 
suspicious straight section visible in the quantum Fred- 
holm determinant (A) in fig. |^ turns out to indeed in- 
dicate a pole, while the Vattay determinant converges 
super-exponentially. 
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FIG. 1. For a analytic hyperbolic map, specifying the con- 
tracting coordinate Wf, at the initial rectangle and the ex- 
panding coordinate z v at the image rectangle defines a unique 
trajectory between the two rectangles. In particular, w v and 
Zh (not shown) are uniquely specified. 



FIG. 2. log 10 |Cn|> the contribution of cycles of topological 
length n to the cycle expansion ^ C n z n for 3-disk repeller. 
Shown are: (o) l/£o, (V) the Gutzwiller- Voros zeta function, 
(□) 1/CoCi) and (A) the quantum Fredholm determinant. Ex- 
ponential falloff implies that l/(o and the Gutzwiller- Voros 
zeta function have the same leading pole, cancelled in the 
1/CoCi product. For comparison, (O) the classical Fredholm 
determinant coefficients are plotted as well; cycle expansions 

for both Fredholm determinants appear to follow the asymp- 

_ 3 / 2 

totic estimate C n ~ A n . A\ symmetric subspace, with 
center spacing - disk radius ratio R : a = 3 : 1, evaluated 
at the lowest resonance, wave number k = 7.8727 — 0.3847 i, 
maximal cycle length n = 8. 



FIG. 3. Same as fig. |, but with R : a = 6 : 1. This 
illustrates possible pitfalls of numerical tests of asymptotics; 
the quantum Fredholm determinant appears to have the same 
pole as the quantum 1/CoCii hut as we have no estimate on 
the size of pre-asymptotic oscillations in cycle expansions, it 
is difficult to draw reliable conclusions from such numerics. 



FIG. 4. Leading resonances in the 3-disk repeller Ai sub- 
space, (a) for the classical Fredholm determinant, and (b) the 
952 leading resonances of the quantum Fredholm determinant 
F qm . Ratio R : a = 6 : 0, cycle expansions truncated at cycle 
length n = 8. 
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FIG. 5. "Difference" between the cfassicaf Fq and F\ 
spectrum for the 3-disk repeller. All resonances so <a 
of Fq that fall within jRe(s ,a) — Re(si,^)| < 0.1, 
|Im(so,a) — Im(si a)\ < 0.08, of a resonance si t p of J\, are 
deleted from fig. m Note that while the non-leading fami- 
lies of resonances so, a almost coincide with the leading si,^ 
resonances, they are not degenerate. With the resonances as- 
sociated with the Fi spectrum removed, families that connect 
nonleading resonances at Re(s) = with the leading part of 
the resonance spectrum for larger Re(s) are clearly visible. 
R : a = 6 : 1, Ai subspace, maximal cycle length n — 8. 

FIG. 6. Complex s plane countor plot comparison of (a) 
the Gutzwiller-Voros zeta function \og\Z qrn (s)\ with (b) the 
quantum Fredholm determinant log |F' 9m (s)j. The border 
of the convergence of the Gutzwiller-Voros zeta function 
agrees with the location of the abscissa of absolute conver- 
gence, given by the J\/2 leading eigenvalue at Re(s) — 0, 
Im(s) = -0.699110157151 . . .. The quantum Fredholm deter- 
minant can be continued at least a factor 3 further down in 
the complex plane. 3-disk repeller, R : a = 6 : 1, A\ subspace, 
maximal cycle length n = 8. 



FIG. 12. The "quantum Fredholm determinant" zeros for 
the baker's map (M5j) computed from all cycles to cycle length 
nmax = 22. A few leading zeros are true zeros; table |n| lists 
the 7 true zeros, the remaining 17 forming a "wedge" in the 
complex plane which slowly drifts toward higher Im(fc) with 
increased cycle truncation length. 

FIG. 13. Same parameter values 3-disk system as fig. 0: (o) 
the quantum Fredholm determinant compared with (□) the 
Vattay determinant. While the quantum Fredholm determi- 
nant is expected to have a pole at Im(fc) = —1.276625955 . . ., 
the Vattay determinant should be entire, and exhibits numer- 
ically faster than exponential convergence. 



FIG. 7. Contour plot of log 10 |-FXs)| for (a) the classical 
Fredholm determinant, and for (b) the quantum Fredholm 
determinant for the Henon map, a = in ([ll|), prime cycles 
up to topological length 12. 



FIG. 8. Contour plot of log 10 |-F(s)| for (a) the classical 
Fredholm determinant, and for (b) the quantum Fredholm 
determinant for the model flow (Wlh with a = 1/2, prime 
cycles up to topological length 12. 



FIG. 9. Contour plot of log 10 |-F(s)[ for the quantum Fred- 
holm determinant; same as fig. ^, but for a wider range of 
Re{k). 



FIG. 10. Same as fig. |8| for the Gutzwiller-Voros zeta func- 
tion. The single leading pole of the Gutzwiller-Voros zeta 
function manifests itself as the border of the convergence; 
only two leading eigenvalues from fig. ^(b) remain within the 
domain of convergence. 



FIG. 11. log 10 \C n \, the contribution of cycles of topolog- 
ical length n to the cycle expansion ^C„z", for the (□) 
"Gutzwiller-Voros zeta function" and (o) "quantum Fredholm 
determinant" for the baker's map (EEI). 
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TABLE I. All periodic orbits up to 6 bounces for the 
Hamiltonian Henon mapping (^) with a = 6. Listed are 
the topological length of the cycle, its expanding eigenvalue 
A p , the variation in the period of the cycle, and its binary 
code. 



rip 


A 




code 


1 


U. / lolb I 0z43oX 1U 


-O.DU /DZOZlOOlU / 





1 


-U.zyozo4bozoy X 1U 


U. 2 (42918851 1(4 


1 


2 


-0.98989794855 x 10 


/~\ O O O O ■"> <~> O <~1 O <~1 O <~1 ■ ") 


10 


Q 

o 


-U.lolyU / 2 / oy / X IU 


-U.zUDUllozyoooo 


1 nn 
1UU 


3 


0.55896964996 x 10 


0.5393446629166 


110 


4 


-0. 10443010730 x 10 


-0.8164965809277 


1000 


4 


0.57799826989 x 10 


0.0000000000000 


1100 


4 


-0.10368832509xl0 3 


0.8164965809277 


1110 


5 


-0. 76065343718 xlO 4 


-1.4260322065792 


10000 


5 


0.44455240007x 10 4 


-0.6066540777738 


11000 


5 


0.77020248597xl0 3 


0.1513755016405 


10100 


5 


-0.71068835616xl0 3 


0.2484632276044 


11100 


5 


-0. 58949885284 xlO 3 


0.8706954728949 


11010 


5 


0.39099424812 xlO 3 


1.0954854155465 


11110 


6 


-0.54574527060 xl0 b 


-2.0341342556665 


100000 


6 


0.32222060985 xlO 5 


-1.2152504370215 


110000 


6 


0.51376165109xl0 4 


-0.4506624359329 


101000 


6 


-0.47846146631 xlO 4 


-0.3660254037844 


111000 


6 


-0.63939998436 xlO 4 


0.3333333333333 


110100 


6 


-0.63939998436 xlO 4 


0.3333333333333 


101100 


6 


0.39019387269xl0 4 


0.5485837703548 


111100 


6 


0.10949094597xl0 4 


1.1514633582661 


111010 


6 


-0.10433841694 xlO 4 


1.3660254037844 


111110 



TABLE II. The "quantum" eigenvalues obtained from the 
Gutzwiller-Voros zeta function Z qrn (s) and the quantum Fred- 
holm determinant F qm {s) for the baker's map. Only the lead- 
ing eigenvalues are expected to coincide. The digits listed cor- 
respond to those unchanged between the cycle length n — 23 
and n — 24 truncations. 





Zqm{s) 






k 


Real part 


Real part 


Im. part 


1 


0.34051779186 


+0.3405177918632516 




2 


-0.39 


-0.3500689652275 




3 




-0.7175716624 




4-5 




-1.099684 


± 0.115118 


6-7 




-1.59 


± 0.227 



18 



